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Abstract 

We continue to study Lax (L-A, U-V) pairs (LP) joint covariance with respect to 
Darboux transformations (DT) as a classification principle. The scheme is based on a 
comparison of general expressions for the transformed coefficients of LP and its Frechet 
derivative. We use the compact expressions of the DT via some version of nonabelian 
Bell polynomials. It is shown that one more version of Bell polynomials, so-called 
binary ones, form a convenient basis for the invariant subspaces specification. Some 
non-autonomous generalization of KdV and Boussinesq equations are discussed in the 
context. Trying to pick up restrictions at minimal operator level we consider Zakharov- 
Shabat - like problem . The subclasses that allow a DT symmetry (covariance at the 
LP level ) are considered from a point of view of dressing chain equations. The case 
of classic DT and binary combinations of elementary DT are considered with possible 
reduction constraints of Mikhailov type (generated by an automorphism). Examples of 
Liuville - von Neumann equation for density matrix are referred as illustrations . 



1 Introduction 

If a pair of linear problems is simultaneously covariant with respect to a DT, it generates Back- 
lund transformations of the corresponding compatibility condition. In the context of such 
integrability the joint covariance principle JT| may be considered as a classification scheme 
origin. In this paper we examine realizations of such scheme looking for possible covariant 
form and appropriate basis with a simplest transformation properties. It is important to note 
that the proof of the covariance theorems for the linear operators incoporates the so-called 
(generalized) Miura transformation having the form of (generalized) Riccati equation. We 
give and examine here the explicit form of the equality in both general and stationary case, 



look as well. It gives additional nonlinear equation that is automatically solved by DT 
theory and used for generation of t-chain equations |§ . We show how the form of the covari- 
ant operator may be found while some kind of Frechet derivatives of the operator coefficients 
and the transforms are compared. 

Further, it is clear that the choice of DT type defines the class of covariant operators. 
The "classic" DT generalization for the polynomials of a differentiation operator results in 
transformation formulas first obtained in Matveev papers , the form of DT we use below 
is from 0. It is also possible to do it in an abstract way: i.e. to define elementary DT (eDT) 
for any projector (idempotent) element in a differential ring or module 0. The abstract 
differentiation is naturally included in the scheme |]|. The appropriate Zakharov - Shabat 
(ZS) problem which contains some elements as "potentials" and "wave functions" (WF) is 
studied. The results open the way for DT theory applications in arbitrary matrix dimension 
(including infinite - dimensional or strictly operator case). It is shown || that the sequence 
of n such transformations innxn matrix case with usual differentiation by a parameter give 
the standard dressing formulas of ZS problem . 

In refs. ||, || eDT for direct and conjugate problems were used to obtain the binary DT 
(bDT) for matrices and in || - for three projectors with applications to N-wave interaction, a 
step from j5| upto the general case |IU| . The symmetric form of the resulting expressions for 
potentials and wave functions (WF) make almost obvious the heredity of reduction restric- 
tions H and underlying automorphisms ([□]] of a generic ZS problems. In JT2J an application 
to some operator problem (Liuville - von Neumann equation) is studied. 

In the Sec. 2 we examine the general non-Abelian version of the operator transformations 
constructed by means of the above-mentioned Bell polynomials generalization. The resulting 
transformation formulas do not contain double summation that is convenient for the further 
analysis as in the context of classification as on the way of dressing chain reformulations |Tj| , 
]14jl , |15| | and solution generation. In the last part of the Sec. 2 we list few such polynomi- 



als. We show then the equation for the intermediate (sigma) function that naturally appear 
during the Darboux theorem proof. This equation generalizes Riccati equation (Miura trans- 
formation) for the celebrated KdV theory for general dynamics - it may be called generalized 
Burgers equation. 

In the next Sec 3 we concentrate our attention on the algorithmic derivation of Lax pairs 
starting from the idea of covariance. We use more simple nonabelian case to show main 
features of the scheme. The point we demonstrate is based on a comparison of the formal 
Taylor expansion of the transformed operator (by means of Frechet derivative notion) of a 
Lax pair and its formal DT expression. The comparison gives differential equations which 
solutions produce the explicit expressions for the Lax pair coefficients. 

Sec 4 introduces binary Bell polynomials and study the transformation of them via the 
differential forms incorporated. We also prove invariance of some combinations of such BBP 
with respect to DT. 

Next section starts from covariance of generalized ZS equation. Two versions of dressing 
chain are constructed for stationary spectral problems with nonabelian spectral parameter. 
We are trying to pick up relations that do not depend on matrix representation [|17[] . Simple 
examples are studied by closures of the chain equations on the operator level. 

The section 6 is devoted to the nonabelian Lax pairs from the point of view of bDT covari- 
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ance. The Abelian version of binary transform (commuting "matrix elements") reproduces 
results of the dressing method based on the matrix Riemann-Hilbert problem [0. 

Here we mainly treat the dressing chain equations, sending readers to [18], where the 
equation 

-iX t =[X,h{X)\, (1.1) 

h(X) - analytical function, is studied. Further generalizations for essentially nonabelian 
functions (e.g. h(X) = XA + AX, [A,x] ^ 0, Hl2| ) are considered in |Tj|, where abundant 
set of integrable equations is listed. The list is in a partial correspondence with |2TJ, and 
give direct link to solutions via the dressing chains. The papers also contain examples of 
"self-scattering solutions" with discussions of possible applications [EU . 



2 Darboux transformations in terms of Bell polynomi- 
als. 

The solution of the problem of a polynomial linear differential operator left and right division 
produces a version of the Darboux theorem. The theorem was reformulated for the classical 
DT 

ip[l} = Dip-atf;, (2.1) 

in terms of some generalization of the Bell polynomials || that give compact expressions for 
transformed coefficients of the operator. We reproduce here the results with links to further 
applications, namely with equations for the element "a" that appear inside the theorem 
proof. 

Let us start from the linear operator 

N 

L=J2 a nD n , 
n=0 

and the evolution equation (flow) 

& = 14- (2-2) 

Here the operator D may be a differentiation by some variable and ip t is the derivative with 
respect to another one (see || for generalizations). The transformation of the solutions of the 
equation is taken in the standard form (2.1), where a = (_D</>)</> -1 with a different solution of 
(2.2) incorporated <pt = L<p. We would now promote the convenient formulation of Matveev 
theorem [[| 

Theorem. The coefficients of the resulting operator 

N 

L[l] = J2 a n [l]D n 

n=0 

are defined by 

ajv[l] = un, 
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and for all other n, by 

N 

a n [l]=a n + [akBk,k-n + ((Dak) - crak)B k -. 1)k -x-n] (2.3) 

k=n+l 

that yields a covariance principle. 

It means that the function ^[1] is a solution of the equation 

V t [i] = L[iMi], 

where L[l] has the same structure and order as L. The proof of this statement incorporates 
the equation that links o and coefficients of the operator L. It may be derived inside a 
factorization theory || and will be useful further 

N 

a t = Dr+[a,r], r = ]T a k B k , (2.4) 

k=i 

where B k are nonabelian Bell polynomials (see, e.g. ||22|| ). In the stationary case one has 

°t = M~ l - M^M' 1 = 0- (2.5) 

The functions B m ^ n are introduced in [|]. 

We reproduce here the definition and some statements about them. 
Definition 

B n , (<T) = l, n = 0,1,2,..., 

and recurrence relations 



B njk (a) = B n _ lik (a) + DB n „ l)k _ x (<7), k = 1, n - 1, n = 2,3,.... (2.6) 

B n>n (a) = DB n _i )n -i(o) + B n (a), n = 1, 2, . . . . 

define the generalized Bell polynomials B m n . 
The following formula is extracted, 

5 n , n _ fe+ i(cT) = J2 [I) B n,n-i(cr) D*~V, fc = T7n, n = 0, 1, 2, . . . ; 

i=k \ / 



S n+1 ( < 7)=X;S n , n _i( < 7)DV, n = 0,l,2,... (2.7). 
it gives the link between standard (nonabelian) Bell polynomials and the generalized ones: 

n 

B n+1 (a)=J2B n , l (v)D n - i a, n = 0,1,2,.... 
Evaluation of the first three generalized Bell polynomials by the definition gives 
B n> x(a) = a; B n , 2 (ff) = o 2 + n Do; B n , 3 (a) = o 3 + no' o + {n - 1) o Do + (™) D 2 o] 
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3 Towards the classification scheme. Joint covariance 
of L-A pairs 

The basis of the formalism we introduce is elaborated starting from [|IJ and the compact 
formulas with the generalized differential (Bell) polynomials from the previous section. Note 
again, that it is valid for nonabelian entries as well, the coefficients a n , solutions of the 
equation (1) <p and ip we consider as matrices or operators. Let us however start from the 
scalar case. 

To prepare the explicit expressions for such work and show details we would setup a couple 
of examples of the theory. Let us deliver a very simple analysis for better understanding of 
the sense of the integrability notion we introduce. First of all we notice that the "elder" 
coefficients, with n = N and n = N-l are transformed almost trivially. It follows that in 
general the functions rather do not play the role of potentials or unknown functions for a 
nonlinear equation of a compatibility condition. 

If N=2, the general transforms (2.3) reduces to 

a 2 [l] = a 2 = a(x,t), 

ai[l] = ai(x, t) + Da(x, t) 

ao[l] = ao + Dai(x, t) + 2a(x, i)Da + o~Da(x, t) (3.1) 

Keep in mind that we touch hear only the scalar (abelian, more precise) case. One can see 
that the explicit form of the transformations really shows a difference between the coefficients 
a(x, t), a,i(x, t) that transform without solutions account and the ao = u(x,t) to be an un- 
known function of a forthcoming nonlinear equation that we call potential in the context of 
Lax (L-A-pair) representation. One may easily recognize the KdV case here. Namely when 
a = const, a\ = , ao plays the role of the only unknown function of the KdV equation. So 
we may formulate an 

observation: The abelian case N = 2 is the first nontrivial example of covariant operators 
set with coefficients a 12 that depends only on x and additional parameter (say t) but its 
transforms contain the only functions a\ y2 , hence, to be called trivial. The transformation 
(DT) for u is given by the last equation of (3.1) and depends on both ai j2 and solutions of 
the equation (2.2) via a. 

By the next order example we would show more details especially when pairs of operators 
are analysed simultaneously. Especially it is important when both operators are functions of 
the only potential. Let us take N = 3. 

Mi] = 63, 

b 2 [l] = b 2 + Db' 3 

b^l] = 61 + Db 2 + 3b 3 Do + aDb 3 

Ml] = 6 + Db x + aDb 2 + (a 2 + (2Da))Db 3 + 3b 3 (aDa + D 2 a) (3.2) 

Consider (3.1) and (3.2) as coefficients of a Lax pair operators both depending on the 
only variable u simultaneously and suppose the coefficients of the operator are analytical 
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functions of u and its derivatives with respect to x. Now we change t — > y and L — > L\ in 
the equation (2.2) corresponding to the case (3.1) and left the parameter t (L — > L 2 ) for the 
case (2.5), forming a Lax pair 

ip y = Lxi). (3.3) 
*l>t = L^. (3.4) 
To produce the KdV case generalization we go to a stationary in y solution of (3.3). 

(j)y = \(f), 

where A is a constant. 

Let us also recall the KdV case. Then the stationary version of (2.4) for N = 2 is 

2 


that reads as 

a 2 + a' + u = c = const. (3.5) 

Note that the equation (2.5) for N = 3 is still valid for this a = (f)^" 1 , if is a solution to 
the Lax pair system (3.3,3.4). If in the equations of (3.2) we restrict ourselves by the case of 
b 2 = and b 3 = b = const, we arrive at the second equation of the KdV Lax pair. 

Next, returning to the general case and taking into account the triviality of b 3 = b(x,t) 
and b 2 in the above-mentioned sense, the first non-trivial potential is 

&! = F(u,u',...). (3.6) 

Suppose further that the covariance principle is valid, or, in shorten arguments notations, we 
address to an equation for this function F. 

6i[l] = F(u[l]) = F(u + Da x + 2aDa + aDa) = F(u) + Db 2 + 3bDa + aDb. (3.7) 

Analyticity of F gives a possibility to use Taylor series expansion in the left side of (3.7) 

F(u[l\) = F(u) + F u (2aDa + Da x + aDa) + F Du (...) + .... (3.8) 

In fact we compare the transform (3.7) with the Frechet differential (3.8) of the F. The 
equation holds identically if the coefficients by the a, Da and the free term are the same. 
Introducing F u = c(x, t) yields 

2ac = 36, (3.9) 



or 



with additional conditions 



*W = (3.10) 



cDa x = Db; (3.11) 
cDa = Db (3.12) 
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Plugging c from (3.9) into (3.12), one go either to 3D(lna) = 2D(lnb) and, integrating, 
obtains 

b = a 3 / 2 Cl (t), (3.13) 

or to Da = Db = 0. In the last case (3.12) is valid with arbitrary c, or independent b(t) and 
a(t). While (3.11) yields the equation for a\ for both cases 

3D ai == 2aDb/3b (3.14) 

with arbitrary Ci(t). 

Next conditions follow from the last equation of the system (3.2), i.e., if one introduces 
next analytical function G and denote 

bo = G{u,u',...), (3.15) 

the transformed bo gives 

G(u + Da l + 2aDa+aDa) = G(u)+G u (Da 1 + 2aDa+aDa)+G Du D(Da 1 + 2aDa+aDa) + .... 

(3.16) 

The DT transformation formula for the potential u is obviously used. The DT for the last 
coefficient 6 , see (3.2), yields 

6o[l] = G(u) + Db x + aDb 2 + (a 2 + 2(Da))Db + 3bD(a 2 /2 + Da). (3.17) 

We should account now the general version of the Miura transformation (3.5) that has the 
form 

2 

a n B n — u + ai<j + a(a 2 + Da) = /i, 

o 

by which we would express the a 2 in (3.17). Doing this and equalizing the expressions (3.16) 
and (3.17) yields 

D — - + aDb 2 + ((fi-u- a x a)/a + (Da))Db + 3bD[(fi -u- a x o)j2a + (Da)} = 

G u {D ai + 2aDa + aDa + G Du D(D ai + 2aDa + aDa) (3.18) 
The equation (3.18) gives for the coefficients: by D 2 a 

G { Du)2a = 3b, (3.19) 

by Da, taking (3.19) into account 

G u 2a + 9b(Da)/2a = Db - 3b ai /2a, (3.20) 

by a 



3b 

G u Da + —D 2 a{x, t) = Db 2 - aja - 3bD(a 1 /2a), (3.21) 
2a 



and the free term is 



D ^T + ((A* - u)/a)Db + 3bD[(n - u)/2a\ = G u D ai + 3b{D 2 ai )/2a. (3.22) 



From (3.19) and (3.20) we express G u : 

G u = Db/2a - 3ba 1 /Aa 2 - 9b(Da) /4a 2 . (3.23) 
If G u is nonzero, from (3.21) follows 

3b 

(Db/2a - 3ba x /Aa 2 - 9b{Da) / Aa 2 )Da + — D 2 a = Db 2 - a x /a - 3bD(a 1 /2a) 

The free term (3.22) gives 

Db . . 3bDa .Db 36ai 9bDa.^ , . „, , ^ 9 . ,„ ,„„., 
u- + n(Db/a -_ = (_-_- _ )Dai{x ,t) + 36(^0/2- (3.24) 

When u is linear independent from <r and derivatives, and we do not account higher terms in 
thr Frechet differential, the only choice Db = kills the term with u, and (3.24) simplifies 

D 2 a x - a 1 (Da 1 )/2a = 0. 

The condition Da=0 as the corollary of (3.12) is taken into account. The equation (3.21) 
also simplifies 

Db 2 - ax/a - 3b(Da 1 )/2a = 

and integration gives the expression for b 2 - 

Another possibility is G u = 0. It gives 9b(Da)/2a = Db — 3ba\j2a instead of (3.23). The 
free term transforms as 

U ~2~^ + ^ Dh / a ~ 3bDa / 2a2 ) = +3b{D 2 ai )/2a. (3.25) 

and by the same reasons gives the conditions Db = Da = 0. It further means a± = and, 
finally, from (3.21), 

Db 2 = 0. 

Hence this case contain KdV equation with the possible a(t) , b(t) , b 2 (t) . 

Remark 1. The results for one isolated equation (3.3) contain rather wide class of 
coefficients (in a comparison with the joint covariance of (3.3) and (3.4). Namely, the a, a\ 
are arbitrary functions of x,t. It may be useful for construction of potentials and solutions 
(e.g. - special functions) for the linear Schrodinger and evolution equations of the one- 
dimentional quantum mechanics [Q. The KdV case may be described separately (denote 
further Df = f): 

G u o' + G u ,a" = 36(1 - a)?/ /4a 2 + 3ba" /Aa 
The only choice is possible, if we consider a, a', u' as independent variables 

G u = 0, 

G u > = 3b/ Aa 

or, taking into account the condition of zero coefficient by u', a = 1, one arrives at 

G(u,u',...) = 3bu'/A. 
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The result leads directly to one of equivalent Lax pairs for the KdV equation. 

Now we would consider the next natural example with interchanged equations defining a 
spectral problem and evolution. Hence we should start from the equation (3.3) for the third 
order spectral operator, using (3.4) as evolution with N=2. 

We would restrict ourselves by the case of 6 2 = and the autonomous b 3 — 1 , b± — 
const = b, b = u. The DT yields 

6i[l] = 6i + Sba' (3.26) 

Ml] = b + b[ + 3b(a x a + a xx ). (3.27) 

As it was shown by analysis of the third order operator, see (3.21-25), the covariant 
spectral problem has the form 

if>xxx + Wx + Gip = \ip (3.28) 

The second (evolution) equation of the case is : 

ipt = -ix - wip (3.29) 

If one consider (3.3) ( specified in (3.26) and (3.27)) or (3.28,29) as the Lax pair equations 
both with coefficients depending on the only variable w. Suppose again that the coefficients 
of the operators are analytical functions of w and its derivatives (or integrals) with respect 
to x. If one wants to save the form of the standard DT for the variable w (potential) the 
analysis just similar to the previous example give for the constant a 

b = 3w/2 + a. (3.30) 

Compare with the formula (3.25), or = w + 2o~ x . Then the transformation for the 
potential w follows from the last equation of the system (3.27), i.e., 

G[l] = G + 3w x /2 + 3{a 2 /2 + a x ) x . (3.31) 

we see that further analysis is necessary due to the possible constraint (reduction) existence. 
Then, again similar to the previous case of KdV the covariant equations (3.27,28) are accom- 
panied by the following (Burgers) equation. 

a t = -(a 2 + a x ) x - w x (3.32) 

for the problem (3.28) and 

a 3 + 3a x a + a xx + bo + u = const, (3.33) 

see (2.5), compare with (3.5) that was Riccati equation (stationary version) for the second 
order spectral problem corresponding to KdV. If one would use the equation (3.32) in (3.31), 
the time-derivative of w appear. Moreover, the further analysis shows that the case we study 
need to widen the functional dependence in u, namely we should include not only derivatives 
of w with respect to x, but integrals (inverse derivatives) as arguments of the potential. So, 
for the equation (3.28) let us introduce next analytical function G and denote 

G = G(d~ l w, w, w x , ...d~ l w u w t , w tx , ...), (3.34) 
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the first terms of Taylor series for (3.34) read 

G(w + 2a x ) = G(w) + G Wx 2a xx + G d -i Wt 2a t + (3.35) 

where we show only terms of further importance. The DT transform (24) after substitution 
of (3.31) gives 

G(w) + 3w x /2 + 3(-<7 t - w x )/2 + 3a xx /2. (3.36) 
Equalizing (3.36) and (3.35) one obtains 

G Wx = 3/4; G 9 -i Wt = -3/4, 

That leads to the exact form of the Lax pair for the Boussinesq equation from [2] for the 
choice of a = — 3/4. 

Remark 2. As one could see we cut the Frechet differential formulas on the level that 
is necessary for the minimal flows. The account of higher terms lead to higher flows (higher 
KdV, for example). 

Remark 3 In the second example we obviously restrict the description by the autonomous 
case within the special choice of DT and do not consider the class of gauge equivalent equa- 
tions. The general case with account of higher derivatives will be presented elsewere in 
forthcoming publications. 

We would finish with the 

Theorem. The Darboux covariant Lax pair (3.3,4) w ^h the stationary equation (3.3) 
where the operator L\ is given by (3.1) has the following coefficients: 

a 2 = a(t), 

ai = ai(x,t), 

is a solution of the equation Dai — a\/Aa = c 2 (t) with arbitrary c 2 , 

ao = u. 

The evolution (3.4), which form is presented by (3.2), contains 

b 3 = b(t), 

b 2 = J aidx/a + 3bai/2a, 

h = 3b/2a, 
b = Sba^/Aa 2 + 3b(Du)/2a. 
The pair produce the generalization of the KdV equation 

Lu = [Li, L 2 ] 

that is solved by standard DT. 
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4 Binary Bell polynomials under Darboux transfor- 
mation 

The binary Bell polynomials (BBP) Y that we are going to use are defined in terms of the 
exponential differential polynomials: 

Y mx , nt (f) = e~ f d™d?e f (4.1) 
by means of the link between ^-polynomials and the standard Hirota expressions f2A 



D p x D q t G' ■ G = (d x - d x ,) p (d t - d t ,) q G'(x, t)G(x', t') . (4.2) 

x'=x,t'=t 

The connection is given by the identity: 

v — w — v — W V 4- w w — V 

y mx> ntM = «p[— ] exp[— 2— «p[^-] • exp[— ] (4.3) 

They inherit the easily recognizable partitional structure of the Bell polynomials (BP): 

y x (v) = v x , y 2x (v, w) = w 2x + vi, y x , t {v, w) = w xt + v x v t , 

y-zx{v, w) = v 3x + 3v x w 2x + vl, ■■■ (4.4) 

The expression of the BBP contains the structure we studied in previous sections. There- 
fore the covariance of the expressions can be examined in similar way. To check the conjec- 
ture about covariance of a BBP one should confirm the coincidence of the expressions of such 
polynomial generated by (4.3) and transformations of functions (2.3) . Let, as the simplest 
example, consider the second order polynomial from the point of view of preceding analysis 
reproducing the addition formula from [EH]. We get 



y2 X {v, v + q) = P B 2 (v) + P 2x B . 
Here the symbols Y 2x ,P 2x (q) are also defined in [p4j|,p5| and B n are usual BP. 

y2x(v, v + q) = ip^Lip = *l)~ l (il) xx + q 2x ip) 

Thus, as it was shown in the Sec. 3, the polynomial contains D-covariant expression for the 
operator L (N=2) with the potential q 2x . It means that the spectral problem for KdV case 
may be expressed by means of the BBP, or 

Y 2x (v,v + q) = X. (4.5) 

This simple example shows how to find the basis in which it is convenient to make the 
transformations. It turns out easier than the direct test and allows to forward the business 
for arbitrary N. 

The general expression for the BBP, transformed like (2.3) [EJJ is the addition formula 



E(N/2) 



N 



y N x(v,v + q)= ]T )P 2p xY iN -2p) x (v) (4.6), 



p=0 



.2p 
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where E(N/2) is the integer part of N/2. After the logarithmic linearization (16), 

E(N/2) / AT \ N 

y Nx (v,v + q) = Y, (m- )P(N-m)MD m iJ = £ a^D m ij, (4.7) 

m=N-2E(N/2) \ iV m ) m=0 

the expression now coincide by the form with the right-hand side of a Lax equation with the 
operator L. We see the coefficients of the DO in (4.7), which DT is determined by (2.3). The 
extracted coefficients are 

if N — m is even, otherwise = (for odd N — m. Or in other words, are nonzero only 
in the set of m G {iV, N -2,...,N - 2E(n/2). Substituting the coefficients from (4.7) with m 
= N - 2p into (2.3), one obtains 

N ( N \ 

= «m + Y f [( N _ ■) ( l P (N-j)x( ( l)[ B j,j-m ~ <7Bj_ij_i_ ro ] 
j=m+l \v JV J/ 

+P(N- j )x(q)'B j - ld - 1 - m }. (4.8) 

The terms in the sum of (4.8) are nonzero also if N - j = 2p is even, it is marked by the prime 
in the sum. The definition, convenient formulas and explicit expressions to the generalized 
Bell polyinomials B n m are given in the Sec. 2. 

For more illustration of the technique let us take again the minimal polynomials and 
evaluate the r.h.s of (4.7) and (4.8). First take an even N and begin with 

a7v_ 2 [l] = ajv-2 + a N (B Nt2 - aB N - 1:1 ) = a N - 2 + Na'. 
We have taken into account that = Pq = 1 Compare with the addition formula 

P 2x (q + 2v) = (^j p 2x(q) + Q p 2x(2v) = q 2x + 2v 2x = q 2x + 2a'. 

We see that in the case of N = 2 the polynomial is invariant. 
The example of odd polinomial N =3 yields 

y 3x {v, v + q) = y 3x (v) + 3q 2x v x = ^"\D 3 + 3q 2x D)^fj (4.9) 

with the known transform of the operator inside (4.9). There are two possible scopes of an 
incorporation of the polynomials into Lax pairs. Either it belongs to evolution equation or 
results in a spectral problem of third order and leads to SK or Boussinesq equations. 
Let us first comment the KdV case, so the spectral problem 

y2x(v,v + q) = X. (4.10) 

introduces a potential q 2x to be a dependent variable for the KdV equation. We should 
rewrite now the evolution equation in terms of binary Bell polynomials. Starting from (4.10) 
one notice that the corresponding coefficient of (4.10) are defined by the equation (4.10): 

Pi P o (3) 

a 3 =1, a x = 3gy , 
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The transforms of them is determined by (4.9) or one may also use the explicit expressions 
from (4.5). We should write, for example, 

3af[l] = 3af + 3a. 

It follows that the potential (dependent variable) should be choosen as 

flS = ?xx/2 (4.11) 

to fit the transformation law for q xx . If one examine the transformation of the coefficients of 
the general operator (e.g.(4.6)), one notices that the "youngest" coefficient bo is changed by 
necessity. The addition depends on a potential and a. So the zero value of Oq could not be 
generally preserved. We arrive to the 

observation: the isolated y 3x is not generally covariant. 

Therefore there is a stright necessity to incorporate different polynomials to support the 
covariance. So we take the combination 

c 3 y3 X (v,v + g (3) ) + c 2 y 2x (v,v + q {2) ) + Cl y lx (4.12) 

with different potentials q{. If we restrict our choice by the unique potential that enter the 
spectral problem (4.5) and note that the choice of g^ 3 ^ = q/2 is fixed by (4.11), we get the 
link of and q from (4.5) that gives 

q ( V = q + 3c 3 q x /(4c 2 ) (4.13) 

it means that constants q are arbitrary and influence only the form of the Lax pair and the 
compatibility condition. 

Proposition. The combination (4-12) by the condition (4-13) is covariant with respect to 
DT. If c 2 = C3 = 1, the compatibility condition of (4-5) - (4-12) is equivalent to the equation 

u t = 3uu x /2 + u xxx /4 

that is nothing but one of KdV forms. 

Many of other integrable equations may be incorporated into the scheme to be developped. 



5 The example of nonabelian Zakharov-Shabat prob- 
lem 

Here we outline the general scheme of the dressing chain derivation in the nonabelian case. 
Let us start from the evolution 

a ty + a l Dty = ty t (5.1) 

that contain the simplest form of the polynomial L(D)-operator. The case is the nontrivial 
example of a general equation (2.3) with non-abelian entries. We would respect \I/ or other 
(necessary for DT construction) solution $ as operators. The equation enter to the Lax 
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pair of integrable nonlinear equations ( e.g. - Davey-Stewartson) while its stationary versions 
produce Nonlinear Schrodinger case. It is also interesting by itself: quantum one- dimensional 
Pauli or Dirac equations after multiplication by the appropriate matrix takes the form of the 
equation (5.1). 

The potential may be expressed in the terms of a from the stationary version of the 
equation (2.4,5), when a t = 0,. Namely, introducing the iteration index i, we have the link 

Da l + [oq, a*] = -Dai ~ K, v&i- (5.2) 

the connection is linear, but contains the commutator. Let us denote ad a x = [<J,x]. Then 

a l = (D - ad ai )~ l {-Dai - [a u 0^.} (5.3) 

The existence of the inverse operator in (5.3) need some restriction for the expression in {,} 
brackets, the expression should not belong to the kernel of the operator D + ad Si . In the 
subspace, where the Lie product is zero, the equation (5.2) simplifies. The DT is also simple. 

4 +1 = 4 + [°i> *i] ( 5 - 4 ) 

Note that a x is not transformed due to the general formulas (2.3). Substituting the link (5.3) 
for i,i+l into (5.4) one arrives to the chain equations. One also could express matrix elements 
of do in terms of the elements of the matrix a and plug into the Darboux transform (5.4) 
separately. 

Let us give more details of the construction in the stationary case. Let us denote ao = u, 
a\ = J, choosing the case DJ = and $ correspond A,/i. Note that there are two 
possibilities for stationary equations from nonabelian (5.1): either 

= A* 

or 

% = *A (5.5) 

and the first of them leads to the essentially trivial connection between solutions and poten- 
tials from the point of view of DT theory Iffi]] . Namely, if 

J$ x + u$ = 

then 

o = Q^- 1 = J- x (ji-v) (5.6) 

and DT do not contains eigen functions. It means, for example, that starting from constant 
potential one never obtain x-dependent u by iterations. 
In the second case one writes 

a = $ a ,$- 1 = J- 1 ^^- 1 - u) 

and DT takes the following form 

u i+l = u i - J- 1 ^, J] + J-V> J] = J-WJ - J- 1 ^, J], (5.7) 
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where it is denoted s = 1 , here and further iteration number indices omitted. The 

potentials u l may be excluded from the equation (2.4) for this case 

a t = Dr + [r, a], 

with 

r = Ja + u = s (5.8) 
. The stationary case, after the plugging u l from (5.8) and returning indices, gives 

s l+1 = s l + Ja l+1 - a\J (5.9) 

and the link of s with a a is very simple 

Ds l + [s\a 1 } (5.10) 

The formal transformation that leads to the chain equations is similar to the (5.3), namely, 
after substitution 

a { = -ad~ 1 Ds i 

into the equation (5.9). 

Further progress in the development of this programme is connected with the choice of 
the additional algebraic structure over the field we consider. It can be useful for the concrete 
representation of solutions of the equation (5.10). For example, if the elements belong 
to a Lie algebra with structure constants C^p, then, after the choice of the basis e a one 
introduces the expansions 

S = 

and 

o l = re- 
plugging into (5.10) gives the differential equation 

+ = 0- 

If one defines the matrix 

B Pa = C^C v (5.11) 

then, outside of the kernel of B 

4 = - B £&- 

By the definition of the Cartan subalgebra C the correspondent subspace does not contribute 
in the Lie product of (5.10). 

Statement. If, further, J belongs to a module over the Lie algebra, Je a = Jp a ^p o,nd there 
exist an external involutive automorphism such that e+ = — e a , C 1 *^ = C^p, e a J = J a pep. 
Then, the chain equation for the variables takes the form 

where the matrix B is defined by (5.11) and the components e a outside of C. Otherwise 

DC, = 0, 
15 



if e 7 G C. 

The system of differential equations is hence nonlinear as the matrix B depends on £, 



The related scheme is developed in Po 

Remark. The scheme may be generalized for a nonstationary equation (2.5). 
equation (5.2) then have the additional term Da\ from the r.h.s. 
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The 



6 Covariance with respect to Binary DT and Dressing 
Chains 

In this section, as in previous, we do not touch directly the joint covariance with respect to bi- 
nary DT (bDT)and classification problem, sending the reader to the publications fl2|,Hl8|,|20 



an, especially to most general [TjJ, where the nobelium functions and rational dependence 
on parameter are admissible. Those papers are devoted to evolution in quantum statistics: 
the equation of Liuville-von Neumann is considered as compatibility condition of stationary 
ZS problem algebraic reduction and equation that introduce a time variable. 
Hence we continue to study stationary version of the ZS problem (5.1) 

JDtf+u¥ = ztf (6.1) 



and dressing chains of bDT (see the Introduction, pi] , compare with e.g. [|K|, the definition 
of which we use here). The algebraic reduction imply D^> = —p^/,Dp = 0, pisascalar.[ 
After the reduction we use the straight corollary of (6.1), changing the notation as in fl2| , 
J — > A, u — > p: 

pAip + pip = z^<p, (6.2) 

and the conjugate problem 

v\A + XP = ZvX, (6-3) 
Let the projector P have a representation 

P = P(X,VVX, (6.4) 



for a given constant idempotent p [10|. The bDT formula connects the elements (operator 
valued potentials - density matrix in the physical context), 

p[i + l] = p\i] + (fM i -vd[P i ,A] (6.5) 

. We would rewrite the conjugate spectral problems [|T^] in terms of P, by means of the 
definition (6.4) 

(j>- f iA)P = z tt P 
P{p - vA) = z v P 

Each equation links the "potential" p and the projector P that plays the role of the element 
a of previous sections. 

pP = fiAP + z^P 

Pp = vPA + z v P. (6.6) 
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The importance of both of them is obvious from the following corollaries 

P(p - pA)P = (p - pA)P 
P(p-uA)P = P(p-vA). 

Let us also note that 

(y - p)PAP = ( Zfl - z v )P (6.7) 

and 



VZp - pz v 
v — p 

One can check that the DT may be treated as similarity transformation 



p p = ^ ^ ( 6 . 8 ) 
v — p 



p[i + 1] = (l + ^^P,)p\A (l + ^—^-Pi) (6-9). 

Multiplying (6.9) by Pj + i from the left side and by Pi by the right side, changing p[n+l] = 
p n one obtains: 

Pn+lPn+lPn = P n +l(l + ~ — -P„)p„-P„(l + fraCU n - p n V n P n ) 

or, using the links (6.6) for the appropriate n, 

P n+ i(u n+1 A + z^ n+ i)P n = P n+ i(l + — — —P n )(p n AP n + z^ n P n ){\ + fracv n - p n v n P n ). 
Simplifying, one arrives to the chain equation. 

{Vn+l ~ Vn)Pn+\AP n + (z^n+l — Z^ n )P n+1 P n = 0. (6.10) 

or, similarly, after the right multiplication by P n +i and left one by P n , 

(.Pn+1 Pn) Pn-^Pn+1 (Z(j,,n+1 ^/x,n )P n P n+1 = 0. (6.11) 

These equations are recurrences and define the solution via parameters. The dependence 
of the parameters on the additional variable t is introduced by the "master" equation (see 

my- ' 

iP t = {AP - PAP) jp- [PA - PAP) jv. (6.12) 

The "minimal" closure P n +\ = XP n = P„A, is possible only if A = 1 for the P n is a 
projector. In this case 

Pn+l Pn P-i Zv,n+1 ^i>,n %v 

So the equations (6.10,11) are valid authomatically. Other simple closure 

Pfi+l O'nPn &n 

with a n , b n e K leads to P n+1 = P n - (1 - (-l) n /2. 
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Next possibility is when some element A do not commute with P n . For example 

Pn+l Pn-^-t ^fi,n+l ^-n^fi 

that for (6.11) leads to 

where, S n = P n A, is proportional to some projector, or, by recurrence, it means 

P n = Po (t)A n 

where the function po(t) is defined from (6.12). 
Finally, we arrive to a conjecture: 

The chain equations (6.11,12) may help to construct particular solutions of the compati- 
bility condition for the Lax pair (see jl^/j if each P n is a solution of the equation (6.13). 



7 Conclusion 

We conclude that our approach 

a) give a connection between a class of a bilinearizable NLPDE, and its Lax form. The 
class is formed by a certain combination of binary Bell polynomials that satisfy the covariance 
principle. 

b) produce automatically a technique of integration of the equations that belong to this 
class. The underlying Darboux covariance of the both linear operators permits to generate a 
rich set of solutions. Between them are multisolitons, positons, periodic and rational solutions 
with abundant combinations that introduce their interaction. [ 

We add that a platform for a nonabelian generalization is prepared. The joint covariance 
principle application presumably cover the items of classification that appears inside most of 



symmetry approaches [pSf|,pB]| and contain the direct link to solutions via dressing chains - 
iterations of DT. 

Mention also the generalization of the theory of small deformations of iterated transforms 
II with respect to intermediate parameters that appear within the iteration procedure of 
bDT . Infinitesimal invariance generated by extended Backlund transformation was consid- 



ered in |30| . The perturbation formulas allow to define and investigate generators of the 
corresponding group, that is a symmetry group of a given hierarchy associated with the ZS 
problem under consideration. 

We continue to develop the dressing method for linear and soliton problems. The pro- 
cedure of dressing includes both solutions and potentials of underlying linear system. The 
other reductions and generalizations may be introduced if in the ring there exist an involutive 
automorphism. 

There are abundant possibilities for generalizations based on recent results of V.B. Matveev 



[27j, where the DT transformation operations is defined via some automorphism at rings. 

Acknowledgments 
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